Direct current superconducting quantum interferometers with asymmetric shunt 

resistors 
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We have investigated asymmetrically shunted Nb/Al-A10a;/Nb direct current (dc) superconduct- 
ing quantum interference devices (SQUIDs). While keeping the total resistance R identical to a 
comparable symmetric SQUID with R' 1 — R^ 1 + R^ 1 , we shunted only one of the two Josephson 
junctions with R = Simulations predict that the optimum energy resolution e and thus 

also the noise performance of such an asymmetric SQUID can be 3-4 times better than that of its 
symmetric counterpart. Experiments at a temperature of 4.2 K yielded e f» 32 h for an asymmetric 
SQUID with an inductance of 22 pH. For a comparable symmetric device e = 110 h was achieved, 
confirming our simulation results. 

PACS numbers: 85.25.CP, 85.25.Dq, 74.25. F- 74.40.De 



The transport characteristics and noise performance 
of direct current (dc) superconducting quantum interfer- 
ence devices SQUIDs having symmetric Josephson junc- 
tions has been intensively studied from the 1970's. Nu- 
merical simulations of the Langevin equations describ- 
ing the SQUID dynamics reliably helped to understand 
the modulation patterns V($> a ,I) and the low-frequency 
voltage noise power Sv(&a, I), where V is the dc volt- 
age across the SQUID, / is the bias current and $ a is 
the applied flux. With the flux-to- voltage transfer func- 
tion V$ = \dV/d$ a \, one obtains the flux noise power 
5$ = Sv /V£ or energy resolution e = S$/2£, where L 
is the SQUID inductance. For an optimized device one 
obtains in the limit of small thermal fluctuations an en- 
ergy resolution e = (8 — 9)ksTL/R for an inductance 
parameter /3x = 2IqL/<&q somewhat below 1 [HQ. Here, 
Io and R respectively denote the junction critical cur- 
rent and resistance. $ is the flux quantum. Although 
e can be very low - for example, in Ref. a value of 
~ 3 ft has been reported at 4.2 K for a 2pH device - 
one may ask whether or not it still can be improved by 
introducing asymmetries in the junction parameters or 
perhaps by adding new elements to the SQUID. Early 
simulations have shown that asymmetries in the junc- 
tion critical currents and resistances can enhance V$ , al- 
though for the prize of asymmetric U($ a ) patterns [![. 
It has also been predicted that an additional damping 
resistor can enhance V$ 0, H| . Several works addressed 
junction asymmetries and additional damping resistors 
in more detail [sj-Tllj] . with the result that the transfer 
function can be increased and flux noise be decreased. 
The above investigations, however, explored only a very 
limited range of parameters and often addressed devices 
where the symmetric counterpart was far from optimum. 

Let us start with a theoretical analysis, using the stan- 



dard Langevin equations [l[ where the Josephson junc- 
tions are described by the resistively and capacitively 
shunted junction modelJH, [ljj. With i — I/Iq the nor- 
malized currents through the junctions k — 1, 2 are given 
by 



-±j = p c (l±a c )5 k +(l±a r )Sk+(l±ai)($) sm(5 k )+i N ,k 

, a c , ct r and c*i denote the asymmetries in capacitance, 
resistance and critical current respectively. The junction 
critical currents are io.fc = -^o(l i a i)i their resistances 
i?/(l ± a r ) and their capacitances C k = C(l ± a c ). '±' 
refers to junctions 1 and 2, respectively. 8 k denotes the 
phase of junction k, j = J/Iq is the normalized circulat- 
ing current in the SQUID loop and (3 C — 2-kIqR 2 C/$o is 
the Stewart-McCumber parameter. Dots denote deriva- 
tive with respect to normalized time r = Qq/2ttIoR. The 
normalized noise current ijy,& has a spectral power den- 
sity AT, with T = 2nk B T/Io^o. The S k are related by 



5 2 -5i = 2tt$/$ + ttPlU 



(2) 
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where <E> is the total flux through the SQUID. L = 
L\ + L2, where L\ and L2 are the inductances of the 
two SQUID arms, related to the inductance asymmetry 
a L via L k = L(l ± a L )/2. 

From Eqs. ([IJ and ([2]) one obtains the normalized 
dc voltage v = V/IoR, and thus the current voltage 
characteristic (IVC) by taking the time average of u = 
(5i + 52)/2. From a Fourier transform of u one obtains 
the normalized correlation functions s v = Sv2ttIqR/^, 
s$ = S*I R/(2$>ok B T) and e = s /2r/3 L . 

The quantity we are interested in most is the optimized 
normalized energy resolution e op t, where optimization is 
done for some or even all SQUID parameters. Recently, 
we have performed a systematic optimization of the noise 
performance of the rf SQUID, optimizing all of its pa- 
rameters [Til . [lFjj . We now apply the same procedure 
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to fully optimize e of the dc SQUID, with respect to 
i, (f> a — ^q/^oj PliPc, oti, a c and a r , i.e., for a given value 
of one or some of these parameters we find all others so 
that e op t is minimized. The inductance asymmetry oil 
does not appear in the above list, since, for a given bias 
current, it only causes a phase shift in v(<j> a ) and also the 
noise correlation functions. Since in practice the junction 
capacitance is always nonzero, the McCumber parameter 
/3 C should be as large as possible to obtain large values of 
IqR. It turns out that /3 C values below 0.8 are uncritical 
in the sense that the other parameters can be tuned so 
that e op t attains its minimum value irrespective of /3 C . 




FIG. 1. Optimized normalized energy resolution e op t vs. (a) 
resistance asymmetry a r and (b) junction critical current and 
capacitance asymmetry a% — a c . Fixed parameters are F = 

0. 01 and fi c = 0.7. For all data points Pl has been varied. 

Fig. QTa) shows e op t vs. a r for on — a c — (full 
circles) and for variable a, = a c (open circles). Fixed 
parameters are j3 c = 0.7 and T — 0.01. The parameters 

1, 4> a and (3l have been varied to minimize e. We have 
used on — a c , having in mind junctions where both Iq 
and C scale with the junction area and R can be chosen 
independently by shunting. The value for T was chosen 
as typical for operation at 4.2 K. For on = 0, e op t vs. a r 
decreases from ~ 1.6 to 0.7 for a r —> 1. In case of variable 
on the minimum e opt is about 0.4, i.e., a factor of 4 lower 
than the energy resolution of a comparable symmetric 
SQUID. Note that each point in the graph corresponds 
to different values of i, (f> a , Pl- We do not list the value 
of these parameters explicitly but note that in all cases 
(3l was in the range 0.4 — 0.5. on = a c ss 0.3 — 0.5 was 
found in case these parameters were varied. 

In Fig. []Jb) we show e op t vs. on = a c for variable Pl 
and fixed (i c = 0.7 and T = 0.01 and a r = 0.99. The low- 
est values of e op t are achieved for on near 0.5. For lower 



values of on, e op t monotonically increases. In particu- 
lar e op t is achieved when on and a r have the same sign, 
i.e., the junction having the lower resistance should have 
the higher Iq. We have obtained similar results, giving 
almost the same lowest values for e op t, also for higher 
values of T (up to 0.1). e opt is thus a robust quantity. 

In dimensioned units e = e • 2^>oksT/ IqR. To maxi- 
mize Iq for /3l ~ 0.5, L should be as small as possible. 
Then, to maximize R and keep j3 c below 1, C should be 
as small as possible, which for a given capacitance per 
area means to keep the junction area as small. If junc- 
tion asymmetries are considered, given a constant critical 
current density, the size of the weaker junction is presum- 
ably limited by the fabrication process and, to obtain an 
on of, e.g., 0.3, the average junction area is increased by 
about 40% from its minimum value. This basically com- 
pensates the gain in e op t- Asymmetries in on are thus 
not necessarily helpful. Thus, below we discuss an ex- 
perimental design having on = a c = 0. 

We also note that for a r very close to 1, e op t increases 
slightly again with increasing a r for a, = a c = 0. This 
is related to chaotic dynamics which appears in some 
ranges of i and <fi a . Below, we will address this issue in 
comparison to experimental data. Nonetheless, if small 
values of e op t can be retained for a r — > 1, the easiest way 
to realize the corresponding SQUID experimentally is to 
"move" the shunt from junction 1 to junction 2, leaving 
junction 1 unshunted and junction 2 shunted with a resis- 
tance R/2. In the following we discuss the performance 
of such a Nb/Al-AlO^/Nb SQUID and compare it to sim- 
ulations, as well as to the performance of a corresponding 
symmetric SQUID. 



asymmetric 



symmetric 




FIG. 2. (Color online) Optical image of the (a) asymmet- 
ric and (b) symmetric SQUID. The Josephson junctions are 
labeled by "JJ" and the shunt resistors by R and R/2. (c) 
Readout scheme of the asymmetric SQUID using a SQUID 
amplifier operated in flux-locked loop. 

Figs. (2fa) and (b) show optical images of an asymmet- 
ric and a symmetric SQUID. The SQUIDs have been fab- 
ricated using a Nb/Al-AlO^/Nb technology based on op- 
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tical photolithography. By a combination of reactive ion 
etching (RIE) employing CF 4 and 2 and ion beam etch- 
ing (IBE) the ground-electrode was defined. The subse- 
quent definition of the junction area was done by RIE and 
anodic oxidation in an aqueous solution of (NFLjBsOs 
and C2H6O2. Before the definition of the vias (RIE - 
IBE) the connecting bridges for the anodization between 
the individual SQUIDs were removed. The following def- 
initions of the resistor, insulation layer and wiring layer 
were all done using a lift-off technique. For the resistor 
material a 76 nm thick Palladium layer was deposited, re- 
sulting in a sheet resistance of 1 fi/sq at T = 4.2K. The 
~ 300 nm thick SiO insulation layer was deposited using 
thermal evaporation while the samples were mounted on 
a water cooled copper plate (T proccss < 26 °C). After in- 
situ pre-cleaning, the final Nb wiring layer connecting 
the junctions, vias and shunt-resistors was dc-magnetron 
sputtered at room temperature. 

Transport and noise measurements were performed at 
T = 4.2 K in a magnetically and electrically shielded en- 
vironment. Dc characteristics (IVC, V($ a ), critical cur- 
rent 7 c ($ a )) were measured in a standard four-point con- 
figuration, using low noise current sources and a high 
impedance room temperature voltage amplifier (RTA). 
For the noise measurements the RTA was not sensitive 
enough. Thus, V was preamplified with a commercial 
SQUID amplifier [l6[ having a 60pV/Hz 1 ^ 2 resolution, 
operated in a flux-locked loop with ac flux bias at mod- 
ulation frequency / mo d = 256 kHz. The SQUIDs were 
operated open loop at fixed I and $ a . V was mea- 
sured by connecting the input coil of the SQUID am- 
plifier in parallel to the SQUID. A 5 f2 resistor Ri was 
in series to the input coil, as shown in Fig.[5Jc). Due to 
the low input impedance of the amplifier, the current 7 CS 
from the current source divides into / and the current 
Ii through the input coil. At given 7 CS , I varies when 
changing $ a , affecting V$ and thus the determination of 
the energy resolution. Using Kirchhoff's laws and the 
condition for flux-locked loop operation, IiM\ = //Mf, 
with Mi (Mf) being the mutual inductance between the 
amplifier SQUID and the input (feedback) coil, one ob- 
tains Ii = (M{ /Mi) (Vf/Rf). Since Mi, M f and R f are 
constants, Ii can be determined by measuring Vf. To de- 
termine U($ a ) for constant I a software control-loop was 
implemented adjusting J cs such that, for each value of $ a , 
I = I cs -Ii was fixed. The ratio (M;/M f ) slightly differed 
for different samples and was adjusted for each device 
until the bias corrected V($ a ) curve measured with the 
SQUID amplifier fitted the corresponding V(<fr a ) curve 
measured with the RTA; all other measurements for a 
given device were then performed with fixed (Mf/Mi). 

Fig. [3] (a) shows IVCs of the asymmetric SQUID. Solid 
black line is for $ = 0, solid gray line for $ a = 0.5 $o- 
The critical current is I c « 2I = hi + I02 = 62.0 ^A 
and for R/2 we obtain 0.57 f2, yielding I C R = 35.3 ^V. 
One notes that, in contrast to IVCs of symmetric de- 
vices, the IVC of the asymmetric SQUID exhibits several 
structures, including regions of negative differential resis- 




FIG. 3. (Color online) Dc characteristics of the asymmetric 
SQUID: (a) IVCs at $ a = (solid black line) and at <Ev = 
0.5 $0 (solid gray line). Lower right inset: 7 c ($ a ). Upper 
left inset: IVC at $ a = 0.843 $0. The arrow indicates the 
voltage for the lowest energy resolution. Theoretical curves 
are shown by dashed lines, (b) V(& a ) (solid black line), for 
/ = — 76.4 /j,A ... 75.9 fiA (in 4.9 [iA steps). Corresponding 
theoretical curves are shown by dashed lines. 



tance. These structures are reproduced in simulations, cf. 
dashed lines. The negative differential resistance in fact 
separates a high-current regime having chaotic dynamics 
from a more stable low-current regime. For the simu- 
lation we have used parameters 0l = 0.675, f3 c — 0.27, 
T = 0.0065, a r = 0.999, a, = a c — 0. These parameters 
have been inferred partly by fitting the IVC, but also 
by fitting J C ($ Q ). The corresponding data for J c ($ a ) are 
shown by solid black lines in the lower right inset of Fig. [3] 
(a). The dashed line in this graph shows the calculated 
curve. Finally, from (3l and Jo we obtain L = 21.7 pH 
which is close to the design value of 23.9 pH. The upper 
left inset of Fig. [3] (a) shows by solid black line the IVC 
taken at $ a = 0.843 $o- For this particular flux value the 
best energy resolution was found at the voltage indicated 
by the arrow. The dashed line is a calculated curve, using 
the parameters given above. A family of curves V($ a ) 
for variable I is shown in Fig. [3Jb). Experimental data, 
for different I are shown by solid black lines. The corre- 
sponding calculated curves, shown by dashed lines, fit the 
data reasonably well, showing that the dc characteristics 
of our device can be understood by the SQUID Langevin 
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equations. In the V($ Q ) curves one notes that the slope 
dV/d& a is very steep for $ a > 0.5 $ ) hi fact reaching 
maximum values of about 1.2mV/$o near / = 56 fiA. 
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FIG. 4. (Color online) Electric transport and noise vs. $ a for 
the asymmetric SQUID at optimum / =56 fiA, measured with 
SQUID amplifier (dots) in comparison with numerical simula- 
tions (line plus symbol): (a) Voltage across the SQUID; solid 
black line shows corresponding curve measured with a high 
impedance room temperature amplifier, (b) Transfer func- 
tion V<s> —dV/d^a- (c) Voltage noise Sy 2 , experimental data 
averaged between 100 Hz < / < 3kHz (white noise regime), 
(d) Absolute (left) and normalized (right) energy resolution. 
In (c) and (d) the noise of the SQUID amplifier has been sub- 
tracted; within the shaded area the SQUID amplifier noise 
was above the noise of the asymmetric SQUID, resulting in 
large errors when calculating e and e. The vertical dotted line 
indicates the position of minimum e. 

Our central results are shown in Fig. |4] Figure EJa) 
shows U($ a ) for the optimum bias current of 56 /i A. 
Dots represent the experimental data, as taken by the 
SQUID amplifier. For comparison the solid black line 
represents the corresponding data taken by the RTA. The 
two curves are well on top of each other, justifying our 
method of correcting the bias current in measurements 
using the SQUID amplifier. The line with symbols is the 



theoretical curve obtained from numerical simulations, 
which agrees well with the experimental curves. Figure 
HJb) and (c) , respectively, show by dots the experimental 

1/2 

V$ and the white voltage noise S v , and in comparison 
the corresponding calculated curves (lines plus symbols). 

1 /2 

For experimental data, S v has been averaged between 
100 Hz < / < 3 kHz. The voltage noise of the SQUID 
amplifier has been subtracted. Figure EJd) displays e 
and e, calculated from the graphs shown in (b) and (c). 
For this device the optimum energy resolution is 32 h; in 
normalized units, e opt = 0.52. Surprisingly, the theoreti- 
cal value for e opt = 0.85 is higher. The reason for this is 
an instability in the calculations, appearing as a chaotic 
switching between two nearby voltage states [17]. This 
seems to be absent in the experimental device. The min- 
imum rms flux noise was 133 n^o/Hz 1 / 2 which is also 
quite low. For comparison, for the symmetric SQUID 
having parameters I R = 37.07 ^V, /3 L = 0.74, (3 C = 0.18 
and T = 0.00526 we obtained e opt = 110 ft (e opt = 1.79) 

and Sl /2 = 361n$ /Hz 1/2 , i.e., a factor of 3.4 higher 
e op t than for our asymmetric device. Further note that 
for the asymmetric SQUID f3 c was only 0.27, allowing in 
principle to increase R by a factor of ~ 1.5, potentially 
decreasing e to ~ 20 h. 

In conclusion we have shown that asymmetries in the 
junction parameters of a dc SQUID can help to signifi- 
cantly improve its noise performance over a SQUID with 
symmetric junctions. The most practical way to achieve 
this is to leave one junction unshunted while the other 
junction is shunted by half of the resistance of the shunts 
of a symmetric device. Experimentally, we have found 
a factor of 3.4 improvement of the SQUID energy res- 
olution over a symmetric device with comparable pa- 
rameters. At least in simulations the drawback of the 
asymmetric shunt is to introduce chaotic behavior of the 
SQUID for certain regimes of bias current and applied 
flux. Our experimental device seems to be less sensitive 
to chaos. Thus, the asymmetrically shunted SQUID may 
be useful for applications where ultra-low values of en- 
ergy resolutions are desired. 
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